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Abstract 

To explain the accelerated expansion of the late universe, the 1/R correction to Einstein gravity 
is usually considered, where R is the Ricci scalar. This correction term, if stable, is generally 
believed to be negligible during inflation. However, if the 1/R term is inflaton-dependent, it will 
dramatically change the story of inflation. The entropy perturbation will naturally appear and drive 
the evolution of curvature perturbation outside the Hubble horizon. In a large class of models, 
the entropy perturbation can be made nearly scale- invariant. In Einstein gravity the single-field 
inflation with a quartic potential has been ruled out by recent observations, but it revives when 
the 1/R term is turned on. The evolution of non-Gaussianities on large scale are also studied and 
applied to inflation with 1/R correction. In some specific models, a large non-Gaussianity can be 
naturally generated outside the horizon. Recent study ruled out almost all f(R) models during 
matter dominated phase. Taking this into consideration, we are left with a limited class of model 
which recovers the Einstein gravity soon after reheating. 
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I. INTRODUCTION 



Despite great achievements of Einstein's gravity theory, numerous versions of its modifi- 
cation or extension have been proposed in the last and this century. Some proposals came 
and went, while others were tightly constrained by observations One of the modern mo- 
tivations for modifying Einstein gravity is attempting to explain the accelerated expansion 
of the late universe [If, 0, S]- Rather than introducing a cosmological constant or an un- 
known dark energy, one can explain the cosmic acceleration by designing a modified theory 
of gravity, see typical models in [El, @, 0] for instance. Among the nonlinear modifications 
[1, E|, namely f{R) gravity theories, the most disputed one is a model with 1/R correction 
to the Einstein action 0, 



S 



M 2 



d £y— g \R — — J + / d x^f^gLu-, 



R 



(1) 



though such a term looks bizarre from the viewpoint of effective field theory. 1 

In reference j6[, it was assumed that the 1/R term was negligibly small in the early 
universe but gradually reveals itself as the universe becomes more and more flat (the Ricci 
scalar R gets smaller and smaller) at the late time. However, in general the parameter 
fi may depend on some matter fields and therefore evolve along with the matter fields. 
Indeed, the general coupling between gra vity and the matter sector was considered in recent 
investigations [HI, [H, [3, 14, 15,0, 17, IS]. Especially, the parameter fi can be a function 
of the inflaton field, which will induce a correction term to single-field inflation, 



S 



d x\ 



-fj 



1 -M 2 R+ 9 M 
2 p 2R 



-g^d a( pdp<p - V(<p) 



(2) 



We expect the correction term was not negligible during inflation but decayed soon after 
inflation (during reheating). By fine-tuning the coupling gUp), one may also expect action 
([2]) reproduces ([1]) in the late universe. Please refer to [l9j for a delicate model unifying 
inflaton, dark matter and dark energy with a single field. 

In fact, the action (jSJ) is just a special case of the f(ip,R) generalized gravity, see 0, 



2Q|, |21|, [22|, [23 



2J, 125J and references therein. So we can employ the formalism recently 



developed in [10(] to deal with this model. In section [TT] we will collect the main results 
of [Io| . in a way as general as possible. The non-Gaussianity in f(ip,R) theory has not 
been studied previously and deserves a separate investigation. But a semi- quantitative 
analysis of this problem will be presented in section [TTTl Then these formulas will be applied 
to model (TS]) in section IIV1 where we also find out the conditions for generating nearly 
scale-invariant power spectra. Based on these results, we will in sections [V] and |VT] study 
models with specific potentials, i.e., V(ip) = ^m 2 ip 2 and V(<p) = \tp 4 respectively. The five- 
year Wilkinson Microwave Anisotropy Probe (WMAP5) has ruled out single-field inflation 
with V((f) = \ip 4 in Einstein gravity, because in this model the tensor-to-scalar ratio 
is too large and the the power spectrum of curvature perturbation P n is over-tilted Q]. 2 
Interestingly, our results will show that, in V(cp) = Xip 4 model, ty can be depressed by an 



1 Throughout this paper, we will mainly following the conventions and notations of [10(. Some conventions 
are gathered in the nex t se ction. 

2 See, however, reference [26| for a counter example. 
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order of magnitude while Pji can be less tiled due to the 1/R term, hence the model can 
pass the WMAP5 test. In the past few years, by considering the coupling to matter in high 



redshift 
models II 




it was found that there are instabilities in some branches of f(R) gravity 
Therefore, in section IVI!} we analyze the stability problem for our models 
and its implication to post-inflation evolutions. We will conclude in the last section after 
a few remarks on the possible loopholes and the resulting uncertainty of our calculations. 
In appendix [A] we will derive some useful formulas for three-point correlation functions of 
entropy and curvature perturbations. The formulas developed in section HTT1 and appendix lAl 
are very general and can be applied to other inflation models with weakly coupled multiple 
fields. 



II. INFLATION IN GENERALIZED GRAVITY 



There has been a lot of investigations on perturbation theory in generalized gravity, 
where the Ricci scalar and a scalar field are non-minimally coupled via an arbitrary function 
W 



211, 122, 1231, 124, 123 



S 



d x\ 



-g 



lf( (p ,R)-lg^d a i P d p <f-V(<p) 



(3) 



But most of them are restricted to special cases with only one degree of freedom, although it 
was believed that there should be two degrees of freedom in general 27, 28, 29(1 - In a recent 
research [10|], such an f(ip,R) theory was reanalyzed by incorporating the other degree of 
freedom and the entropy perturbation. Being interested in its implication to inflation, here 
we will gather the general relevant results. In an independent work [3(| , starting with more 
general kinetic terms and more scalar fields, the evolution of the "perturbed expansion rate" 
was calculated for generalized gravity theory using the techniques invented by [3l| . We will 
mainly follow the conventions and notations in lOj. For instance, the signature of metric is 
( — h ++), and we will take 



M: 



8nG, R = Q{2H 2 + H), F 



E 



2HF + F 



(4) 



All of the general results have appeared in [lOj, partly mixing with some special models. 
Nevertheless, it is still helpful to put them orderly in this section. 
First of all, let us define the slow-roll parameters: 



m 
5 3 
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(5) 



Be careful with the notation and sign difference between the slow-roll parameters here and 
those in most literature. The slow-roll conditions are met if the absolute values of these 
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parameters are much smaller than unity. Under the slow-roll conditions, the background 
equations can be approximately written as 



V--f + 3H 2 F ~ 0, 
ip 2 + 2HF — HF ~ 0, 
3H<p - + V^0. 

These equations and the slow-roll conditions also result in the following useful relations 

3Fcp 2 - 4Fip<p 1 



(6) 



E = - 



fI 



E 



2F2 



6- 



2 - ei - -6i, 



6 4 ~ 2771 - ~8 X , 



5 35i5 3 - SxT)! - 4771772 
2 3d 1 — 4?7! 



(7) 



In the longitudinal gauge, the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 
with scalar type perturbations is given by 



ds' 



-(1 + 2<p)dt 2 + a 2 (l - 2ip)S ij dx i dx j . 



(8) 



For generalized gravity, usually <p 7^ ty, so we will have two degrees of freedom after eliminat- 
ing the inflaton fluctuation Sip. The perturbed Einstein equations will give us two coupled 
second-order differential equations for (0, ip). So we say there are two dynamical degrees of 
freedom. But this pair of variables can be traded for (TZ, S) or (v-n, t> 5 ) or (un, Us) by the 
following relations: 



K 
S 

u-r. 



1 

2 

F l~Y 
2F 



2HF + F 
2F(p 2 + 3F 2 
F(2HF + F) 
2F(p 2 + 3F 2 



(2#F + F) 2 
AFip 2 + QF 2 



+ V0 + 



2HF + F 
-^F~ { 



-4) 



a\ 2F(2Fip 2 + 3F 2 ) 



2HF + F 



"ft, ^5 



a\ 2F(2F(p 2 + 3F 2 ) 



AFip 2 + 6F 2 



(0 + VO, 



us 



2HF + F 
asj F(2Fy 2 + 3F 2 ) 
V3(2HF + F) 



S. 



(9) 



We have chosen the normalization for S so that Vn* — fts* when perturbations cross the 
Hubble horizon, as will be given by equation (|14[) . In reference [13], 1Z is interpreted as the 
curvature perturbation while S is interpreted as the entropy perturbation, v-n and v$ are 
the corresponding canonical variables. The interpretation of un and us is less clear, but are 
defined for our convenience, and one may think un as as something akin to the canonical 
momenta (not exactly). In terms of them, the evolution equations of perturbations take the 
form 



u 



'Ilk 

+ k 2 u Sk + m 2 s a 2 u Sk + ak 2 u nk = 0, 



+ k 2 u nk + m 2 n a 2 u nk + (3u Sk = 0, 

7 2 ,2 2 ,7 2 r. 



Sk 



(10) 
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whose coefficients 



m 



H 2 



2 
TZ 



2e 1 -r] 1 , (3 ~ sign(^)aifv^i - 2ei, 



a ~ sign((^)-ai?v5i - 2e i; 

As is well known, couple is trouble. This also applies to the coupled equations (jTUl) . 
To make our analysis simple, for perturbations inside the Hubble horizon, we will always 
disregard the coupling terms controlled by (3 and a (decoupled approximation). This enables 
us to get a rough estimation but also induces some uncertainties. Imposing an appropriate 
quantized initial condition at k ^> aH, we find an analytical solution under the decoupled 
approximation, 

«K* = *e«*-teV^HW(-kT)<h», with v\ = \ - ^, 
4^2 4 H 2 

^ = -^-t)f^«(-k)e a , with^ 2 2 = i-^|. (12) 
2V6A; 4 ii z 



Here {e^, e^fc} is the orthonormal basis 

{e a k,e/3k>) = 8 a p8(k - k'), a,P = TZ,S. (13) 

The normalization of Fourier modes TZ k and S k is exhibited by equation flAlj) in appendix 
[A1 If m^/H 2 ~ and m 2 s /H 2 ~ —2, then the power spectra at the horizon- crossing = ail 
are nearly scale-invariant, 

k 3 |2 # 4 



2tt 2 ' 4vr 2 



2tt 2 ' 4tt 2 5 



(14) 



The condition m^/H 2 ~ is trivial due to the slow- roll conditions. But m 2 s /H 2 ~ —2 
puts a nontrivial constraint on viable models. Throughout this paper, an asterisk means the 
quantities take their horizon-crossing value. Since we have neglected the coupling between 
curvature perturbation and entropy perturbation inside the horizon, their cross-correlation 
is negligible at the Hubble-crossing, 

k 3 

Vc* = ^(n k *,s k *)~o. (15) 

When crossing the horizon, the spectral indices are 

rin* - 1 = n s * - 1 = 4ei* - 2r] U . (16) 

Unlike the single-field inflation in Einstein gravity, the entropy perturbation and curvature 
perturbation are not conserved even well outside the horizon k -C aH. It is more convenient 
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to follow their evolution in terms of (TZ, S). If m|/if 2 ~ —2, we have Sk/(HSk) ~ 0{t) 
and hence 

1 {mi 9 3 \ 

S k = /2 S HS k , /is = -- ( -^f + 2 + e l - -5 X - 35 2 + -5A , 



• sign(^) 2{5 x -2e 1 ) f n 5 \ 

= fi n HS h , fi n = — - — y I 2?7i — -di — 04 J • (17) 

Taking fig and fin as constants approximately, its analytical solution reads 

S k = S fc *exp Q" /i 5 M^ = 5fe»e^^-^, 

- ft** = / ^ffSfcdt = [e«M - 1] , (18) 

Ju 

in which A" = \n[a en d/ ct(t)] stands for the e-folding number from time t to the end of inflation. 
As a result, on the super-hubble scale, the power spectra are 

V S * V s *e 2 ^ N *- N \ 



Their spectral indices are 



/4S 



Tin — 1 = iic, — 1 = — -, 

ns-1 = ns*-l- 2fi s , 

fi s \2e^ N *~ N ^-l] , 
n c~ 1 = "** ~ 1 ~ U("-"> - 1 • (20) 

We have defined the entropy-to-curvature ratio in 13 

- - £ < 21 > 

The tensor ty pe perturbation is conserved outside the horizon. Its power spectrum is 
relatively simple [241 ] 

■>ll- 

V T ^-^ = P 5 *(85i - 1660, (22) 

with a spectral index 

2H F 

n ^Jp-HF =2 ^-^ (23) 

Here we have used a normalization different from 10|, |24( to accommodate to the WMAP5 
convention Q]. 
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The above results hold generally for slow- roll inflation in generalized f(ip, R) gravity with 
F > 0, as long as the entropy perturbation is non- vanishing and nearly scale- invariant. For 
details of derivation and explanations, one can refer to [1 Of ] . 

It proves helpful to utilize also the entropy-curvature correlation angle A 
tensor-to-scalar ratio 

Vc = Vt 



32 and the 



cos A 



(24) 



as well general transfer matrix 



K 
S 



1 T-jis 
T ss 



S 



(25) 



In fact, cos A is nothing else but the correlation coefficient introduced in 33] 



III. NON-GAUSSIANITY 



The primordial non-Gaussianity has attracted a lot of attention during recent years. 3 To 
judge whether we can get some interesting large non-Gaussianity before expanding actions to 
the third order and calculating the three-point correlation functions, we can make some semi- 
quantitative estimates using the results of single-field inflation in Einstein gravity theory. 

We start with the calculation of three-point correlation of curvature perturbations 1Z, by 
virtue of (125L 

(^(k 1 )^(k 2 )^(k 3 )> 

= ((^(kO, + T 7 ^«S(k 1 ) # )(^(k 2 ), + 7WS(k 2 ),)(ft(k 3 )* + 7W>(k 3 )*)> 

= ( l S(k 1 ),5(k 2 ), l S(k 3 )*)r 3 iS 

+ [(7e(k 3 ) < ,5(k 1 )»5(k 2 )») + (^(k 2 )»5(k 3 ),5(k 1 ),) + (^(k 1 ),5(k 2 ),5(k 3 ),)] T^ s 
+ [(7e(k 2 )^(k 3 )*<S(k 1 )*) + (^(k 1 )^(k 3 ),5(k 2 ),) + (7e(k 1 )^(k 2 )* < S(k 3 )*}] Tns 
+ (^(k 1 )^(k 2 )^(k 3 ),), (26) 

where an asterisk means the quantities are calculated at the time of horizon-crossing k = aH. 
In the above equation, we assumed the linear evolution of TZ^ and outside the horizon. 
Although this assumption is good enough for our semi-quantitative analysis, in a more 
accurate treatment, one should consider the nonlinear effects. There are two sources of 
nonlinear effects outside the horizon: the Sk term neglected in equation (iT7|) : the time 
dependence of /is and fin- 

As we have mentioned, the curvature perturbation and the entropy perturbation are cou- 
pled. But, under our approximation, their coupling inside the horizon will not be taken into 
consideration in the estimation of magnitude. Because all of these quantities are calculated 
at horizon- crossing, we can treat the adiabatic and entropy perturbations independently. 



For a p artial list, see 
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43 



44 



60] and references therein. We also recommend 61( as an up-to-date brief overview 
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Using the single-field relations [34j, [35 



(^(k 1 )^(k 2 )^(k 3 ),) = (2vr) 5 / 2 ^ ( 
(5(k 1 )^(k 2 )^(k 3 ),) = (2vr) 5 / 2 5( 3 ) (^k, 



YJ4 



3 ' 
i '"i 



S\ 2 



10 



/iVL* (^fc*) 



(27) 



The power spectra on the right hand side are given by ( 1T41) and related by ( 12TT) . while 
the nonlinear parameters /jv^* and can be estimated independently by the single-field 
results. Hence the summation involving these two pure three-point correlations can be easily 
written into a compact form, 



(«S(k 1 )* < S(k 2 )*<S(k 3 )*)Ti s + (^(k 1 )^(k 2 )^(k 3 )*) 



(2tt 



) 5/2 * (3) (5>) (-^) (^nf 1 



(28) 



The contributions of terms like (1ZSS) and (TZTZS) will be a little trouble before we 
know the exact forms of third order action and perform careful calculations. Here we cannot 
determine the value of these terms for generic configuration, but in appendix IA1 there is an 
estimation for the local shape. We find there the three-point correlations involving both 1Z 
and S are proportional to the two-point cross-correlation Vc, whose initial value is negligible 
at the time of horizon-crossing under our approximation. 5 These proportional relations rely 
on the locality of the shape, although it is possible that they could be generalized to other 
shapes by incorporating the k dependence of f]f L and f^ L . Strictly speaking, we have to 
reevaluate their contributions seriously when going beyond the local limit. But, lacking of 
a solid proof, we will still set (7Z*S*S*) = (7Z*7Z*S*) = for all shapes, which will give us 
some satisfactory results in semi-quantitative estimation. In the previous section, we made 
a decoupled approximation for linear perturbations inside the horizon. The assumption here 
is just a nonlinear generalization of that linear one. After this assumption, we get on the 
super-horizon scale, 



(^( kl )^(k 2 )^(k 3 )) = (2nf^ ( ) ( -1) (P£) 2 g4 (f%L* + rlAfS L# ) 



ru 3 



At the same time, the left hand side of (129]) can be converted into 



<ft(kOft(k 2 )ft(k 3 )) = (2tt) 5 /V 3 > PTk, 



.-If* (v n ) 2 

^qJnl {' k ) 



nA 3 ' 



(29) 



(30) 



4 The notation of the momentum modes of the perturbations in [35j are different from here by some different 
choice of normalization in Fourier expansions. See appendix [X\ for details. 



5 However, as argued in [10( along the line of [62] , the cross-correlation may not be negligible, Vc*/Vs* ~ 
0(a) ~ 0(/3), were the coupling terms taken into consideration. Of course, even if we take them into 
account, due to the a or (3 suppression, the dominant contribution is still given by (1Z1Z1Z) and (SSS) 
terms. So our approximation captures the leading order contributions. 
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where the spectrum can be related with the one at horizon- crossing by (|T9l) . as 

V n = V n * + TlsVs, = (1 + r s *Tls) V n *- (31) 

Comparing (1291) . ( l30i) and ( 13T1) . we finally get the nonlinear parameter f^L of curvature 
perturbation on super-Hubble scale, especially at the end of inflation, expressed by some 
parameters at horizon-crossing, 

f _ J NL* + ' S* TZS J NL* /on\ 

Here /jvl* and f NLif are computed at k — aH. In our approximation, the curvature and 
entropy perturbations are evolving independently before that time, so the nonlinear param- 
eters for them at the Hubble-exit can be estimated with the independent single-field results 



36, 37 



5 

Inl* = -^[^* - 1 + /(ki,k 2 ,k 3 )n T ,], 



s 

12f 



NL* 



2£ 

n s * - 1 - — + /(ki, k 2 , k 3 )n T * 



(33) 



where /(k 1; k 2 , k 3 ) is a factor of momentum configuration, with maximum 5/6 in equilateral 
limit and minimum in local limit j36|. £ is a factor related to the normalization of S. 
Corresponding to our normalization (jHJ), it is 



£ = -\—, — =S 3 -r ]l --5 1 . (34) 

In our model, so far we do not know the relation between the entropy perturbation during 
inflation and the one at the matter-radiation decoupling. So there is an ambiguity in the 
normalization of entropy perturbation. In literature of two-field inflation, a convenient 
normalization is usually chosen so that V-r* = Vs* at the Hubble-exit. We follow the same 
normalization. But one should realize that the it is <S/£ rather than S that satisfies the 
simplest form of the consistency relation. Therefore the second consistency relation (|33|) 
takes a relatively more complicated form. Since n-n* — 1 = ns* — 1 and rs* = 1 to the 
leading order, then we have a simplified estimation of (1321) as 



' i^ + Thf + 6#e(l + 7l 5 ) 



TZS 



5(e + 7l5) ;k-l + /(k 1) k 2 ,k 3 W + -j5|^. (35) 



For nowadays observation, the most relevant results are its values in the local and the 
equilateral limits: 



ciocai __ 5(£ + T 7e<g ) _ _| ^Tfts 



-tzs> v- 1 - ' -•-ns 

r >i<->i 5 (£ + Tjzs) ( _ i , 5 \ 5£T4s 

12£(1 + 7^)' + 6 nT *J + 6He(l + T ns 



fz 11 = \ - 1 + + ^ (36) 



All the parameters involved in this formula can be expressed by the slow-roll parameters 
and e-folding number, as in the previous section. We will evaluate the results for specific 
models given below in section [V] and [VT] 
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IV. l/R CORRECTION TO INFLATION 



With the above results at hand, it is straightforward to study the inflation model (T5J), 
where a inflaton-dependent l/R correction is included. The steps are parallel to those in 
loj. Remember in [loT ] a special model with a inflaton-dependent R 2 term was considered. 



But it turned out the inflaton is rolling up its potential in that model. It is a rather tricky 
problem to terminate inflation in "rolling-up" models. So it would be interesting to get a 
"rolling-down" model in f(ip, R) gravity. The model we are going to study has this quality 
to which we will return at the end of this section. 
Comparing ([5} with d3J), it reads directly, 

f(^R)=M 2 R+ 9 M F = M 2 -jL, (37) 

then we obtain the simplified background equations 

V - 3 M 2 H 2 - ^ ~ 0, (38) 
V? 2 + 2HF -HF-0, (39) 

3H ^iip- v «- (40) 



The solutions for equation (138]) are simply 



1 . , V±JV-lgM} 



H = uff^ = ^IM} ' (41) 

We will always take the positive solution (the one with upper sign "+" ) by virtue of the fact 
H 2 > 0. 

We get the following relations: 

p 16p 2 ' 
p, v (p~6M 2 HH = 2Hpe 1 , 

P,^ + P, w 2 - 6M 2 (M + H 2 ), 
3H<p + 3H<p~(?^-V tl ?j 0, 

3H'<fi + QHip + 3Hip ~ (^j^- - V^J (p+^-^-V^ 2 . (42) 
The slow-roll parameters (jSJ) can be expressed in terms of g and V and their derivatives with 
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respect to ip, 



H_^ Pjf M^f 9jf M^ 

H 2 ~ 2p 2 V 8P ,V 



Hep p V 8p 



F 2e lP 3 M 2 / g \ F 1 

<*1 = 777^ TTTT-o ^TT^ [ ~Z > 5 2 = 7717 - e l ~ TT^i, 



FF p, v (16p 2 -^M2) VpV,/ ^ 
F , 2ep(g/p 2 )^ E 3 

HF PA9/P),v HE 2 



F 2epp, 



FF P 2 



vV 



£ e lC2 2 ei M p 2 / g^M 2 

, E 5 3^3 - ^1 - 4t7!?7 2 

* 6 = is " * - 2 s1 + — — ■ (43) 

The "mass squared" for entropy perturbation reduces to 

™>s 5 e 32p 2 / 3 \ 3g, v M*f g \ ( 1 \ 



To calculate non-Gaussianities, we also need 



(gM%\ 6p 2 



(45) 



In order to move on, we should specify the potential V(ip) and the coupling g(<p). This 
is the task of the coming two sections. Here we should mention the condition to make the 
power spectra scale-invariant. There nontrivial condition m 2 s /H 2 ~ —2 is translated now to 
the requirement 32p 2 /(gM 2 ) ~ 6. According to (l4"Tj) . this requirement is easy to satisfy if 
we choose 

4V 2 - M 4 V 2 

y 3M 2 1 J 

We will take this choice in the subsequent sections. The V w term in the numerator is 
necessary, otherwise one would find F = and pin is divergent. Generally we have 4V 2 ^> 
Mp/ 2 , then from (jlOj) it is not hard to get 

3Hp ~ ~V jip . (47) 

In our following specific examples V = \m 2 p 2 or V = A<p 4 (A > 0), so V^p < and the 
inflaton is ro lling down its potential as promised. We also have F > 0, so the formalism 
developed in lCj is applicable here. 
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V. QUADRATIC POTENTIAL: V(ip) = \m 2 p 2 



In this case, the action takes the form 



S 



1 „ m 4 (o? 4 
-M 2 R + — yT - 

2 p 6M 2 R 



M 4 ) 



\!f*d a <pdfN> - \m 2 p 2 



(48) 



When the scalar field p> fades out, this action recovers the gravitational part of action (jl]) 
if (i A = m 4 /3. But as we will see at the end of this section, this is not the case because 
< m 4 /3. 

For later convenience, let us define a new notation 



7 



M 2 
p 

p 2 



This notation is also useful in the next section. In the present case, one can prove 

4M 2 (3M 4 + </) 



3{M 2 + p 2 ) 



(49) 



(50) 



According to this expression, the condition 7 1 is necessary in order to satisfy the slow-roll 
conditions. As will be clear below, this is also the sufficient condition to meet the slow-roll 
conditions. So we can conclude that this model describes the large field inflation. 

All of the slow-roll parameters can be expressed in terms of 7 to the leading order as 



Cl = 
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H 2 ~ 


4 

"3 7 ' 


Vi = 


Hp 7 ' 


Si 


F 8 „ 
HF 3 ' 


5 2 = 


E 

HE ~ 
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~3 7 ' 


S3 = 
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~~ HF ~ 


S 4 


E A 2 
HE 


£2 = 


H 

HH ~ 


8 2 
3 7 ' 


V2 = 


Hp " 47 ' 


Se~- 


E 

— ~ 4 7 
HE 



(51) 



The coefficients f ill I) are 



m 



H 2 



- 5 7, p 




a 



4 /2 7 



ml 



(52) 



When the perturbations cross the horizon, 

3m V* 



r s * = 



256vr 2 M p 6 : 



n n * - 1 = n s * - 1 



16 



-7* 



32 



-7*, 
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_> 
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47* • 



(53) 
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Figure 1: The evolutions of power spectra with respect to e- folding number TV* — N 
after crossing the horizon. This figure is drawn according to the model with action 
(|48|) . From top to bottom: curvature power spectrum V-ji, entropy power spectrum Vs 
(dashed blue line) and cross-correlation power spectrum Vq (dot-dashed purple line), 
tensor power spectrum Vt- All of the power spectra are normalized by Vs*i the entropy 
power spectrum at horizon-crossing. The vertical dotted black lines correspond to 
N* — N = 60. 



If we use the horizon- crossing value to estimate fj,$ and fin outside the horizon, then at 
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Figure 2: The evolutions of correlation coefficient cos A (upper graph), entropy-to- 
curvature ratio r$ (its logarithm, middle graph) and tensor-to-scalar ratio (lower 
graph) with respect to e-folding number N*— N after crossing the horizon. This figure 
is drawn according to the model with action (|48|) . The vertical dotted black lines 
correspond to — N = 60. The horizontal dotted black line corresponds to cos A = 1, 
that is, the totally correlated situation. 
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/•local, equil 
NL 




Figure 3: The evolutions of nonlinear parameters of curvature perturbation with re- 
spect to e-folding number iV* — N after crossing the horizon. The solid blue line 
corresponds to the local limit value f l § c £ l , while the dashed purple line depicts the 
nonlinear parameter of equilateral shape f^f- The vertical dotted black line cor- 
responds to N* — N = 60. This figure is drawn according to the model with action 

(gg]). 



the end of inflation (N = 0), 



Vs* 64 7 * ^ ! 

V S _ -647V, 7 »/9 Vt_ 64 



V s * V s * 3 

V c 3 



*) 



V s * 8\/T* 



-32A r »7»/9 A _ e -32/V»7»/9^ 



16 16 

n n -l = — -7», n s - 1 = —J*, 

nc-l = y7*, n T = --7*. (54) 

To determine the parameter 7*, we make use of the observational constraint on curvature 
spectral index n-n — 1 ~ —0.04, which gives approximately 7* = 3/400. The results are shown 
in figures [T] and figure [2J In figured] we plot the evolution of spectra V-ji, Vs, Vc and Vt, 
which are defined in (|14p . f[T5"j) and (1221) . When drawing the graph, we have normalized them 
by Vs*- Figure [2] depicts the evolution of the correlation coefficient cos A, the logarithm of 
entropy-to-curvature ratio rs and the tensor-to-scalar ratio Tt, defined by (I21"j) and (f24|) . 

At the end of inflation, it can be seen from figure [2] that the entropy perturbation and 
the curvature perturbation are almost totally correlated, and the entropy-to-curvature ratio 
r s is of order 10~ 3 . At first glance, the entropy-to-curvature ratio here can be tested against 
WMAP5 constraint [H as done by jTo|. However, this is a misleading game. What WMAP5 
constrained is the entropy perturbation 

between dark matter and radiation. In our model the entropy perturbation 13 

Scx^-^-0) (56) 
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which is related to the difference between the Newtonian potential and the spatial curva- 
ture if). Firstly, the normalization of S does not match to S. Second, it is unlikely that the 
two degrees of freedom in our model will decay into radiation and dark matter respectively. 
Most probably such an entropy perturbation would seed an anisotropic stress or quadrupole 
moments of photons and neutrinos. Third, the entropy mode may decay after inflation, 
which depends on the detailed mechanism of reheating. Especially, the entropy perturba- 
tion can be erased by thermal equilibrium of matter and radiation before the creation of any 
non-zero conserved quantum number 0, [63|, 64 . 

To estimate the non-Gaussianity, we calculate fl36l) for the present case, 



r local 
J NL 



requil 
J NL 



40 



20 



2048 7 ^ + 135 V2 (l - e 



-32V»7,/9j 3 



64 7 * + 9 (1 - e - 327V *W9) 2 



174087, 3 + 945^ (l - e 



-327V*7„/9j 3 



27 



64 7 , + 9 (1 - e - 32JV *W9) 2 



(57) 



Once again we set 7* 



1 NL 



3/400, then the numerical result gives the nonlinear parameter f\ oc "'' 
and fxjf as functions of the e- folding number A*. Both of them are illustrated in figure 
|3j At the end of inflation, this model will give the nonlinear parameters /jy" 



(local 



11 and 



requil 
J NL 



13. 



Although all of the above predictions (or postdictions) are consistent with observational 
data, this model suffers from a serious problem, as we want to point out here. If we take 
n-ji — 1 = —0.04, then the normalization of curvature power spectrum V-ji ~ 10 -9 re- 
quires m 2 /Mp ~ 10~ 5 . On the other hand, the smallness of cosmological constant requires 
/i 2 /M 2 ~ 10~ 121 in action (j2J). In other words, if one intends to use model (T48|) to explain 
the comic microwave background (CMB) anisotropy, the residual "dark energy" will be too 
large compared with the observed value. For this reason, we conclude the model (148!) with 
a quadratic potential is unattractive. In section IVlIl we will discuss another problem of it. 



VI. QUARTIC POTENTIAL: V(cp) = A<^ 4 
Starting with the action 



S 



-M 2 R + 



4AV(</ - 36M. 



6M 2 i? 



(A > 0) (58) 



the treatment of this model is similar to the previous section, but the result is more encour- 
aging. 

Again we find the necessary condition 7 = M 2 /<£> 2 1 for slow-roll because of the 
relation 



ei + m 



8M 2 (72M p 6 + 96M p V + 16M p V + ^ 
3</? 2 (6M 2 + y? 2 ) 3 



(59) 
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Pr/Ps* 




Figure 4: The evolutions of power spectra with respect to e-folding number TV* — N 
after crossing the horizon. We draw this figure according to the model with action 
(158|) . The upper graph depicts the evolution of curvature power spectrum Viz- The 
middle graph depicts the evolution curves of entropy power spectrum Vs (dashed blue 
line) and cross-correlation power spectrum Vc (dot-dashed purple line). The lower 
corresponds to tensor power spectrum Vt- All of the power spectra are normalized 
by Vs*-, the entropy power spectrum at horizon-crossing. The vertical dotted black 
lines correspond to N* — N = 60. 
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Figure 5: The evolutions of correlation coefficient cos A (upper graph), entropy-to- 
curvature ratio rs (its logarithm, middle graph) and tensor-to-scalar ratio (lower 
graph) with respect to e-folding number iV* — iV after crossing the horizon. This figure 
is drawn according to the model with action (|58|) . The vertical dotted black lines 
correspond to N* — N = 60. The horizontal dotted black line corresponds to cos A = — 1, 
that is, the totally anti-correlated situation. 




To the leading order of 7, we write the slow- roll parameters in the present case 

H 16 (p 8 F QO 2 

E 16 . F 16 . E 16 

62 = he^-y^ 5s = ^"y 7 ' 6i = m^-^ 

H 16 2 '<p 8 E 16 . 

e 2 = — - ~ 7 , 77 2 = — — 7, (J 6 = — - — 7- 60 

HH 3 1 ' ' HCp 3 h b HE 3 ' V ; 
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Figure 6: The black dot is the prediction of quartic model (I58j) for n-jz and tt, where 
we have set 7* = 1/400. It is consistent with the constraint from WMAP5 + BAO 
(baryon acoustic oscillations) + SN (supernovae) [4]. 




Figure 7: The evolutions of nonlinear parameters of curvature perturbation with re- 
spect to e-folding number N*—N after crossing the horizon. The solid blue curve cor- 
responds to the local limit value /^£ a/ , while the dashed purple curve plots the value 
in equilateral limit f^™ 1 . The vertical dotted black line corresponds to N* — N = 60. 
This figure is drawn according to the model with action (|58|) . 



and the coefficients f[TT|) for evolution equations, 



2 



8 2j mo 
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When the perturbations cross the horizon, 



'Pn* — ^Ps* 



5127r 2 M p 6 : 



r s * = 1, n-n* - l — ns*-l — -I67*, 
176 



/As* = — g-7*, Hn* = g V7*, 



^7*- (62) 



If we use the horizon-crossing value to estimate fis and fi-ji outside the horizon, then at 
the end of inflation (N = 0), the power spectra and spectral indices are 

Vs* 484 % ^ j + A ' 

Ps _ -3527V, 7 »/9 _ 256 

' T-> o 7*) 



Ps* V s * 3 

_ 3 c -i76jV t7 ,/9 _ e -176Af,7»/9^ 



7> 5 * 22^ 
n K - 1 = -I67*, 715 - 1 = -g-7«> 

1 128 32 ^ 

From the horizon-crossing to the end of inflation, the power spectrum of curvature pertur- 
bation has increased significantly. In sharp contrast, the entropy perturbation drops down 
exponentially with respect to iV* — N. The cross-correlation between them takes a positive 
value, at first increasing in amplitude and then decreasing. As we have promised, the tensor 
type perturbation is invariant. These results are presented in figure HI 

Now turn to figure [51 Look at the upper graph for the evolution of correlation coefficient 
cos A. Under our approximation, at the time of Hubble-crossing (N* — N = 0), the curva- 
ture perturbation and the entropy perturbation are uncorrelated. But the the subsequent 
evolution makes them almost totally correlated at the end of inflation (iV* — N = 60). 

From the lower graph of figure [5], it is clear that the tensor-to-scalar ratio is depressed 
greatly even at the horizon- crossing time (compared with the Ay? 4 inflation in Einstein 
gravity). Let us take a closer look on this point. Given the choice (|46p . at the time of 
Hubble-crossing, the power spectrum for curvature perturbation and that for tensor type 
perturbation can be written as 

3V 3 V 

Vn * = 32tt 2 m6\/2' Vt * = 2^nq- (64) 

In contrast, the counterparts in Einstein gravity are given by 

V 3 2V 

The difference between ( I64jl and fl65l) explains the smallness of tt at the horizon- crossing 
in our model. Outside the horizon, since the curvature perturbation is increasing while the 
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tensor type perturbation is conserved, the value of tt becomes smaller and smaller. At the 
end of inflation, we have tt — 0.03. This is well inside the constraint of WMAP5 jiJ], as 
illustrated in figure El 

Words are needed here about the entropy-to-curvature ratio r 5 , which is depicted in the 
middle graph of figure Thanks to the exponential decrease of entropy perturbation, the 
value of r$ is of order 10 -4 at the end of inflation, which is smaller than the WMAP5 upper 
bound [H . However, as we have emphasized in the previous section, it is not quite reasonable 
to compare the entropy perturbation here with that in WMAP5 result. A more relevant 
constraint might come from the quadrupole moments of neutrinos. Especially, since the 
entropy perturbation is very small at the end of inflation, we can treat <fi + ip as almost 
constant at that time. From equation (Q, this gives 

^ _ ^ = Fi^HF^Fl 
2Fip 2 + 3F 2 

in our specific model. It is an interesting question to investigate the implication of residual 
difference between spatial curvature and Newtonian potential. But we do not pursue it 
furthermore in this paper. 

The non-Gaussian features can be studied as before. By virtue of the new relations 
between the slow-roll parameters and 7* we have 

220 [212967, 3 + 9^2(1- e - 176 ^>/ 9 ) 3 

rlocal L 

JNL — r „i 2 



3 

110 



484 7 * + 9(1 - e - 176J W 9 ) 2 



5962887? + 207^ (l - e -iraW9) 3 
f e N q f = l — r -rs— (67) 



27 



484 7 * + 9 (1 - e - 176 ^*/9) 2 



One can plot the evolution of /jy£ ' and with respect to iV* — iV in figure At the 

end of inflation, this model will give the nonlinear parameters /^£ a/ — 9 and f^ d — 12. 
Of course, these numbers are just results of semi-quantitative estimation. If we take them 
seriously, we would like to compare them with observational constraints 0]. They perfectly 
satisfy the WMAP5 limit -9 < f^f < 111 and -151 < ffif < 253. 



VII. STABILITY ANALYSIS 

For model ([T]) in the large curvature region, due to the 1/R suppression, the correction 
term has negligible effects in the early universe if \i is small [6]. But recent investigations 



13j, [14| showed that such correction terms may introduce instabilities, hence their effects 
are not negligible even in the high redshift epoch. So it is important to study the stability 
problem 6 in our model Since the inflaton field is evolving, and there is a signature 
change in g(ip) around the end of inflation, we should study this problem during inflation 
and after reheating respectively. 



We are grateful to the referee for putting this problem to our attention. 
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To analyze the stability, there is an indicator B given by formula (17) in [13| and formula 
(2) in [3]. According to the results of 13, 14], the instability resides in the branch of models 
with B < 0. We find for our model of the form (TSl), the indicator 



B 



2g 



dlnR fdlnH 



MIR 2 - g dlna 



dlna 



(68) 



During inflation, for g taking the form (06]) and M 2 V w <C V, we get B ~ 2 > 0, so the 
model is stable. This can also be inferred simply from the fact that g > during inflation. 
As the inflaton rolled down the potential and decayed long after the inflation, for the 



quadratic and quartic potentials, we have V <C M 2 V !tp(p and then 



B 



2V 2 



dlnR fdlnH 



ZR 2 + V 2 W dlna 



dlna 



< 0. 



The inequality is saturated if and only if V jV>v> = 0. This cannot happen when V 
because it always gives V, 



(69) 



12 2 

|m cp 



m > 0. As a result, it does not have a proper matter- 
dominated phase. This is another problem of the model with a quadratic potential, as we 
have promised at the end of section IVII 

But the story is a little different for the case V = Ay? 4 , in which the unstable branch 
with B < can be avoided if the inflaton decayed to the minimum of its potential ip = 
during the reheating era. Therefore, the stability condition for the quartic potential model 
puts a constraint on reheating: the reheating process should to efficient enough to guarantee 
a complete decay of inflaton ip. Such an efficient decay can be realized most easily by the 



instant preheating mechanism [651 ] . After the complete decay of ip during reheating, the 
f((p,R) model of the form (1581) is reduced to the Einstein gravity without any harmful 
instability in subsequent epochs. 

The stability condition makes our models less interesting. Were the stability condition 
ignored, one might drive the acceleration of late universe with the residual non-vanishing 
inflaton, and thus unify the two phases of accelerated expansion of the universe with a single 
non-minimally coupled scalar field. After imposing the stability condition, there is no room 
for such a natural unification. 

Of course, inspired by the so-called mCDTT model in [3], one may replace g(<p) with 
g(ip) + ^M 2 in (T5]), where fi <C M p is a constant independent of <p. Then the newly added 
term will play a role after inflation, exactly recover so-called mCDTT model, which can 
avoid the instability problem. But, since the constant is very small, this term plays no 
role during inflation. Moreover, as was advocated in all f(R) modified gravity during 
the matter phase is grossly inconsistent with cosmological observations. So it seems that 
the only choice for us is to recover Einstein gravity after reheating. 

Therefore, although we started with the f(R) model and the accelerated expansion of 
late universe, due to various difficulties put forward in 11, H, 13, 14], it turns out that our 
model has nothing to do with the f{R) model at the late time. In other words, the survived 
f(p, R) inflation model should reduce to the Einstein gravity after reheating. 



In contrast with [ll|, [l2j , different viewpoints are held by the authors of [66|, |67|, [68 



Interested readers may refine the analysis above by taking [66|, [67|, [68J] into consideration. 



7 We thank S. Odintsov and S. Nojiri for bringing 
informing us about 

0,0- 



G6 



G7 



68| to our attention, and thank C. Corda for 
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According to [6^, the interferometric detection of gravitational waves can provide a 
definitive test for general relativity. In other words, the interferometric detection of gravita- 
tional waves will be a strong endorsement for the modified gravity theories or, alternatively, 
will rule out them. So it would be also necessary to further inspect the /(</?, R) models from 
this angle of view in the future. 



VIII. COMMENTS AND CONCLUSION 

As we have stressed, our analysis throughout this paper is not more than a semi- 
quantitative estimation. Before concluding, we would like to remark on several weaknesses 
and the resulted uncertainties in the above calculations. We can classify them into three cat- 
egories: the decoupled approximation inside the horizon, the linear evolution approximation 
outside the horizon and the slow-roll approximation. 

First, as revealed by equations fflUj) . the curvature perturbation and the entropy pertur- 
bation are coupled inside the horizon. But when writing down the analytical solution ffl2l) . 
we have neglected the coupling terms. As a subsequence, the correlation functions (7^*5*), 
(IZ^IZ^S*) and (7Z*S*S*) vanish only because we have neglected the coupling between 1Z* 
and inside the horizon. All of the power spectra and three point functions at the horizon- 
crossing should receive a correction from the coupling effects. The correction is controlled 
by coupling coefficients a and (3 in evolution equations ffTUl) . This is also a general problem 
for analytical solution of multi-field inflation models. For a more accurate treatment to this 



problem in two-field inflation, please refer to [62 



Second, in deriving the transfer relation (fT9|) . we have neglected the nonlinear effects. As 
mentioned in section 1111} there are two sources of nonlinear effects outside the horizon: the 
Sk term neglected in equation (TT7T) ; the time dependence of us and ji-n. Again, this is also 
a general problem for analytical solution of multi-field inflation models. 

Third, there is an additional source of uncertainty for the model studied in section HVl 
where we have deliberately kept the V w term in the Lagrangian. This is necessary to avoid 
the divergence of power spectrum at the leading order, but it brings some inconsistency 
for our slow- roll approximation. This is clear from equations (151 p and (1601) . in which the 
slow-roll parameters are not of the same order. In principle, this problem should be solved 
by doing the calculations at the sub-leading order in a consistent way. But the background 
dynamics will be rather messy, neither analytical nor numerical method can give it a hand. 

Although the analytic results obtained in this paper are not accurate, it is still meaning- 
ful to take them for rough estimate before painstaking calculation. There are some lessons 
we can learn from it. For the l/i?-corrected inflation, the evolution of entropy pertur- 
bation can dramatically depress the tensor-to-scalar ratio and enhance the magnitude of 
non-Gaussianity. Specifically, if we take the rough estimation seriously, then the single-field 
inflation can be rescued by the 1/R correction, otherwise it would have been excluded by 
observational data. 

The preliminary investigation in [l(| and here raises more questions than answers about 
generalized f((p,R) gravity theories. First, we lack a first principle to write down the 
exact form of f(ip,R) when higher or lower order corrections are considered. Second, all 
of the calculation makes sense only semi-quantitatively, so a more accurate treatment is 
in demand. The formalism we developed is applicable to other cosmological stages and 
scenarios. Especially, it would be interesting to find a unified model similar to but with 
richer phenomena. Third, it is possible that the entropy perturbation in f(ip,R) inflation 
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can seed a tiny quadruple moment of neutrinos, which deserves a detailed analysis. Fourth, 
according to our rough estimate, the non-Gaussianity is large and positive in some models. 
This is observationally interesting and should be studied carefully in the future. 
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Appendix A: THREE-POINT CORRELATIONS OF THE LOCAL FORM 

Before engaging ourselves in calculation, we notice that the definition of IZk and Sk in 
the text is 

/H 3 k r H 3 k 

— -e* kx ^ k (t), 5(x,t) = / — -re^Ct). (Al) 
(2%) 2 J (2n)2 

But in calculating non-Gaussianity, usually a different normalization is followed, 

/d 3 k ~ r H 3 k 

" e*-*fc k (f), S(x,f) = J " e*-*^*). (A2) 

In terms of TZk arid Sk, we have the following relations between two-point correlations and 
power spectra [3J, [38j: 

(^ k i^ k2 ) = —Vnimk, - k 2 ) = ^LVnimk.! - k 2 ), 

(5 kl 5 k2 ) = — VsWOL! ~ k 2 ) = ^LVsimk, - k 2 ), 

(n kl s k2 ) = ^cwsfa - k 2 ) = V^Lcwsfa - k 2 ), 



/rl 3 k 

/ii 3 k 

/rl 3 k 
—V c (k)e^-^. (A3) 

In accordance with the WMAP convention [U, [34j , we parameterize the nonlinearities of 
curvature and entropy perturbations as 

S(x,f) = S L -^ L (Sl-(Sl)). (A4) 
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Here TZl and Sl are linear Gaussian parts of the perturbations. If we take nonlinear pa- 
rameters /jy L and fx L as constants, then this is a local form non-Gaussianity, which can be 
written in the Fourier space as 



TZ(k) = TZ L (k) - -f> NL 



f/ " P 7^(k-p)7t L (p) 



(27T) 



d d xe- tk - x (7ei(x)) 



5(k) = &(k) - 



f/ " P 5 L (k-p)5 L (p) - / rf 3 xe- ik - x (5i(x)) 



(27T) 



(72-1) and («S|) are counter terms to ensure (7£(x, t)) = (<S(x,i)} = 0. 

Using the above relations, it is straightforward to prove equation (1301) and 

(^(k 1 )^(k 2 )<S(k 3 )) 
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C /x 1 rfx 2C /x 3 e- lkl - Xl - ik2 - X2 - ik3X3 (^ i (x 1 )^ i (x 2 ))(^ L (x 1 )5 L (x 3 )) 



6 f w 

NL 



6 f W 
' NL 



c/x!<ix 2 (ix 3 e 



— iki-xi— ik2-X2— ik3X3 



(7e L (x 1 )7e L (x 2 ))(7e L (x 2 )5 L (x 3 )) 
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A 3 



•^(2vr) Jvl— ^| o " ( k i + k 2 + k s) 



A 3 



•^(2vr) JVL— ^| 5 ( k l + k 2 + k 3)- 



7,3 

A 2 



By exchanging 1Z <-> 5, one directly writes down 

(«S(k 1 )5(k 2 )^(k 3 )) 



_ ^( 27r ) /JVL— Ta 7^—^ ( k l + k 2 + k 3 
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/V 3 
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A 3 



_ T7T( 27r ) Jnl^3 P5— d v ; (ki + k 2 + k 3 



1,3 
h 2 



(A5) 



(A6) 



(A7) 



The above derivation is valid when f^ L and f§ L are constants. This is the case for the 
local shape non-Gaussianity. We hope the results can be generalized to other shapes as if 
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fx L and ffj L are k-dependent. But this conjecture is to be proved or disproved by a more 
careful investigation in the future. 
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